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1 Introduction 

This work presents an Eulerian-Lagrangian approach to the Navier-Stokes 
equation. An Eulerian-Lagrangian description of the Euler equations has 
been used in ([[[], @) for local existence results and constraints on blow- 
up. Eulerian coordinates (fixed Euclidean coordinates) are natural for both 
analysis and laboratory experiment. Lagrangian variables have a certain 
theoretical appeal. In this work I present an approach to the Navier-Stokes 
equations that is phrased in unbiased Eulerian coordinates, yet describes 
objects that have Lagrangian significance: particle paths, their dispersion 
and diffusion. The commutator between Lagrangian and Eulerian derivatives 
plays an important role in the Navier-Stokes equations: it contributes a 
singular perturbation to the Euler equations, in addition to the Laplacian. 
The Navier-Stokes equations are shown to be equivalent to the system 

Tv = 2uCVv 



where C are the coefficients of the commutator between Eulerian and La- 
grangian derivatives, and V is the operator of material derivative and viscous 
diffusion. The physical pressure is not explicitly present in this formulation. 
The Eulerian velocity u is related to v in a non-local fashion, and one may 
recover the physical pressure dynamically from the evolution of the gradient 
part of v. When one sets v = the commutator coefficients C do not enter 
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the equation, and then v is a passive rearrangement of its initial value (J7J, 
B) 0? 0)- When ^ the perturbation involves the curvature of the 
particle paths, and the gradients of v: a singular perturbation. Fortunately, 
the coefficients C start from zero, and, as long as they remain small v does 
not grow too much. 



A different but not unrelated approach ( [[ll]] , ( [12j ) is based on a variable 
w that has the same curl as the Eulerian velocity u. The velocity is recov- 
ered then from w by projection on divergence-free functions. The evolution 
equation for w 

Tw + (Vu)* w = 0, 

conserves local helicity and circulation (when v — 0). We will refer to this 
equation informally as "the cotangent equation" because it is the equation 
obeyed by the Eulerian gradient of any scalar <fi that solves Tcf) = 0. The 
variable w is related to v. 

w = (yA)*v 

where A is the "back-to-labels" maps that corresponds when v — to the 
inverse of the Lagrangian path map. A(x,t) is an active vector obeying 

TA = 0, 

A(x, 0) = x. Both v and A have a Lagrangian meaning when v = 0, but the 
dynamical development of w is the product of two processes, the growth of 
the deformation tensor (given by the evolution of VA) and the rearrangement 
of a fixed function, given by the evolution of v. In the presence of viscosity, 
v's evolution is not by rearrangement only. It is therefore useful to study 
separately the growth of VA and the shift of v. 

Recently certain model equations (alpha-models) have been proposed 
([[L3|], |TJ|]) as modifications of the Euler and Navier-Stokes equations. They 
can be obtained in the context described above simply by smoothing u in the 
cotangent equation. Smoothing means that one replaces the linear zero-order 
nonlocal operator u = Pw that relates u to w in the cotangent equation by a 
smoothing operator, u = K a Pw. When v = the models have a Kelvin cir- 
culation theorem. They cannot be models for Eulerian averaged equations, 
the equations describing mean flow in turbulence theory. These Reynolds 
equations are not conservative in the limit of zero viscosity: the fluctuations 
introduce additional stresses, the so-called Reynolds stresses, that preclude 
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conservation of circulation along the mean flow. This is a simple yet funda- 
mental objection to the identification of alpha-models with Reynolds equa- 
tions. The alpha-models might be models of Lagrangian averaged equations, 
an entirely different concept. (Operationally one may think of Eulerian av- 
eraging as long time averages at fixed position, and of Lagrangian averaging 
as long time averages at fixed initial label. However, a fixed initial label has 
no obvious physical meaning when one deals with ensembles of flows.) 

In this paper we consider the Navier-Stokes equations and obtain rigorous 
bounds for the particle paths and for the virtual velocity v. The main bounds 
concern the Lagrangian displacement, its first and second spatial derivatives, 
are obtained under general conditions and require no assumptions. Higher 
derivatives can be bounded also under certain natural quantitative smooth- 
ness assumptions. 

Some of our bounds can be interpreted as a connection to the Richardson 
pair dispersion law, one of the empirical laws of fully developed turbulence 
([|J), that is consistent with the Kolmogorov two-thirds law (||). The pair 
dispersion law states that the separation S between fluid particles obeys 

(\5\*)~ et * 

where e is the rate of dissipation of energy and t is time. This is supposed 
to hold in an inertial range, for times t that are neither too small (when the 
separation is ballistic) nor too large, when viscous and boundary effects are 
important. The law can be guessed by dimensional analysis by requiring the 
answer to depend solely on time and e. Precise laboratory Lagrangian experi- 
ments have recently begun to be capable of addressing Lagrangian quantities 
with preliminary results that seem to be consistent with the Richardson law 
in some ranges (@). If one considers the problem of estimating the pair 
dispersion mathematically one is faced with the difficulty that the predic- 
tion seems to require both non-Lipschitz, Holder continuous velocities and 
Lagrangian particle paths. Our approach allows a rigorous formulation and 
an upper bound 

(\5\ 2 )<3\5 \ 2 + 24:E t 2 + 6e B t 3 , 

that includes a reference to an initial displacement 5q, an initial kinetic energy 
E and a rigorous upper bound on e. The prefactors are probably not optimal. 
The conditions under which such a bound can be obtained are quite general, 
and there are no assumptions. In this paper we chose the case of periodic 
boundary conditions, and body forces that have a characteristic length scale 
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that remains finite as the size of the periodic box is allowed to diverge. The 
bound €b does not depend on the size of box. In many physically realistic 
situations one injects energy at the boundary; in that case one can find e# 
independently of viscosity (0), without any assumptions. 



2 Velocity and displacement 

The Eulerian velocity u(x,t) has three components u l , i = 1,2,3 and is a 
function of three Eulerian space coordinates x and time t. We decompose 
the Eulerian velocity u(x,t): 

,-, , dA m (x,t) , . dnix.t) 

U \x,t) = JL; >- vjx,t) - (i) 

Repeated indices are summed. There are three objects that appear in this 
formula. The first one, A(x,t), has a Lagrangian interpretation. In the 
absence of viscosity, A is the "back-to-labels" map, the inverse of the particle 
trajectory map aHj = X(a,t). The vector 

£(x, t) = A(x, t) - x (2) 

will be called the "Eulerian-Lagrangian displacement vector" , or simply "dis- 
placement" . £ joins the current Eulerian position x to the original Lagrangian 
position a = A(x, t). A(x,t) and £(x,t) have dimensions of length, V A is 
non-dimensional. A pair of points, a = A(x,t), b = A(y,t) situated at time 
t — at distance 6q = \a — b\ become separated by St = \x — y\ at time t. 
From the triangle inequality it follows that 

(St) 2 <3\£(x,t)\ 2 + 3\£(y,t)\ 2 + 3(S ) 2 . (3) 

The displacement can be used in this manner to bound pair dispersion. 

The second object in ([]]), v(x,t), has dimensions of velocity and, in the 
absence of viscosity, is just the initial velocity composed with the back-to- 
labels map (0, ||, 0, pl)| ). We call v the "virtual velocity". Its evolution 
marks the difference between the Euler and Navier-Stokes equations most 
clearly. The third object in (0) is a scalar function n(x, t) that will be referred 
to as "the Eulerian-Lagrangian potential" . It plays a mathematical role akin 
to that played by the physical pressure but has dimensions of length squared 
per time, like the kinematic viscosity. If A(x,t) is known, then there are 
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four functions entering the decomposition of u, three v-s and one n. If the 
velocity is divergence-free 

V-u = 0, 

then there is one relationship between the four unknown functions. 

3 Eulerian-Lagrangian derivatives and com- 
mutators 

When one considers the map x h- > A(x, t) as a change of variables one can 
pull back the Lagrangian differentiation with respect to particle position and 
write it in Eulerian coordinates using the chain rule. Let us call this pull-back 
of Lagrangian derivatives the Eulerian-Lagrangian derivative, 

V A = Q*V E . (4) 

Here 

Q(x,t) = (VA(x,t))-\ (5) 

and the notation Q* refers to the transpose of the matrix Q. The expression 
of Va on components is 

V\ = Qjidj (6) 

where we wrote di for differentiation in the i-th Eulerian Cartesian coordinate 
direction, 

di = V l E . 

The Eulerian spatial derivatives can be expressed in terms of the Eulerian- 
Lagrangian derivatives via 

V* E = An) (7) 

The commutation relations 

[v^,v|] = o, [vi,vy =0 

hold. The commutators between Eulerian-Lagrangian and Eulerian deriva- 
tives do not vanish, in general: 

[Vi,V|] =C mMl VJ. (8) 
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The coefficients C m)k]i are given by 

C m , k , t = {V\(d k £ m )}. (9) 

Note that 

Cm,k;i = Qjidid k A m = V^(V^A m ) = [V^, V^]A m . 

These commutator coefficients play an important role in dynamics. 

4 The evolution of A 

We associate to a given divergence-free velocity u(x, t) the operator 

d t + u-V - vA = T u (u,V). (10) 

We write d t for time derivative. We write T for T u (u, V) when the u we use is 
clear from the context. The coefficient v > is the kinematic viscosity of the 
fluid. When applied to a vector or a matrix, T acts as a diagonal operator, i.e. 
on each component separately. The operator V obeys a maximum principle: 
If a function q solves 

Tq = S 

and the function q has homogeneous Dirichlet or periodic boundary condi- 
tions, then the sup-norm H^Hl^^) satisfies 

||g(-,*)||L°°(dx) < ||g(-,*o)||L°°(<fcO + / \\S(-, s)\\ L oo( dx )ds 

J t 

for any t < t. The operator Y v {u, V) is not a derivation (that means an 
operator that satisfies the product rule); V satisfies a product rule that is 
similar to that of a derivation: 

T(fg) = (Tf)g + f(Tg) - 2v{d k f){d k g). (11) 

We require the back-to-labels map A to obey 

TA = 0. (12) 
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By ( |T^| ) we express therefore the advection and diffusion of A. We will use 
sometimes the equation obeyed by £ 

(d t + u- V-vA)e + u = (13) 

which is obviously equivalent to (0). We will discuss periodic boundary 
conditions 

£(x + Lej,t) = t(x,t), 

where ej is the unit vector in the j-th direction. Some of our inequalities will 
hold also for the physical boundary condition that require £(x, t) = at the 
boundary. 

It is important to note that the initial data for the displacement is zero: 

£(x,0) = 0. (14) 

The matrix VA(x,i) is invertible as long as the evolution is smooth. This 
is obvious when v = because the determinant of this matrix equals 1 for 
all time, but in the viscous case the statement needs proof. We differentiate 
(jjjD in order to obtain the equation obeyed by V 'A 

r(VA) + (VA)(Vu) = 0. (15) 

The product ( V A) (Vm) is matrix product in the order indicated. We consider 



TQ = {Vu)Q + 2vQd k (VA)d k Q. (16) 

It is clear that the solutions of both (|15D and flUf ) are smooth as long as the 
advecting velocity u is sufficiently smooth. It is easy to verify using ( |i~T| ) that 
the matrix Z = (VA)Q — / obeys the equation 

TZ = 2uZd k (VA)d k Q 

with initial datum Z(x,0) = 0. Thus, as long as u is smooth, Z(x,t) = 
and it follows that the solution Q of (|l^) is the inverse of VA 

The commutator coefficients C m ^ enter the important commutation re- 
lation between the Eulerian-Lagrangian label derivative and T: 

[r,vy = ivC m ^ k E v m A (17) 
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The proof of this formula can be found in Appendix B. 

The evolution of the coefficients C m ^i defined in @ can be computed 
using ([15p and (0): 

r (c m>fe;i ) = -($4n)v^(«,)) 

-(d k (ui))C mjl .i + 2uC jMi ■ di {C mMj ) . (18) 
The calculation leading to ([T8|) is presented in Appendix B. 

5 The evolution of v 

We require the virtual velocity to obey 

Tv = 2uCVv + Q*f. (19) 

This equation is, on components 

T v (u, V)vi = 2vC m)k -id k v m + Qjifj. (20) 

The vector / = f(x,t) represents the body forces. The boundary conditions 
are periodic 

v(x + Lej,t) = v(x, t) 
and the initial data are, for instance 

v(x,0) = u (x). (21) 

The reason for requiring the equation (|19|) is 

Proposition 1. Assume that u is given by the expression (|7[j above and 
that the displacement I and the virtual velocity v obey the equations ftH\ ) and 
respectively Then the velocity u satisfies the Navier- Stokes equation 

d t u + u ■ Vii — vAu + Vp = f 

with pressure p determined from the Eulerian-Lagrangian potential by 

Id 2 

T u (u, V)n + — + c = p 
where c is a free constant. 
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Proof. We denote for convenience 

D t = d t + u-V. (22) 

We apply A to the velocity representation ([!]) and use the commutation 
relation 

[D t ,d k ]g = -{Vu)*Vg. (23) 

We obtain 

AK) = m D t A m )) v m + (d t A m ) D t v m - d t fJtf + D t n 
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We substitute the equations for A (|13|) and for v (|19|): 
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A(« l ) = -oii M-L- + D t n + (a,(z/AA m )) t> m + 



(«9^ m ) {z/A Wm + (2^ fc (W)*A^ + /, )} . 

Now we use the facts that 

(d^Q^ = 8 {j 

(Kronecker's delta), and 

d k (V£)* u = d k (VAy u = d^A 1 ) 

to deduce 

D t (u l ) = -di f M- + An) + /• 

+i/(A«9 i A m K + K3^ m )A Wm + 2^(5^')^ 
and so, changing the dummy summation index I to m in the last expression 

AK) = + An) + uAmA m )v m ) + U 

Using (|lD we obtained 
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A(« l ) = ^A^ - J— ! z/An + A^ + /, 



and that concludes the proof. 



Observation The incompressibility of velocity has not yet been used. This is 
why no restriction on the potential n(x, t) was needed. The incompressibility 

V • u = (24) 

can be imposed in two ways. The first approach is static: one considers the 
ansatz (P and one requires that n maintains the incompressibility at each 
instance of time. This results in the equation 

An = V • (VA)*v) . (25) 

In this way n is computed from A in a time independent manner and the 
basic formula (P can be understood as 

u = P((VA)*v) (26) 

where P is the Leray-Hodge projector on divergence-free functions. The sec- 
ond approach is dynamic: one computes the physical Navier-Stokes pressure 

p = R i R j (u i u i ) + c (27) 

where c is a free constant and Ri = (— A)~5<9j is the Riesz transform for 
periodic boundary conditions. The formula for p follows by taking the di- 



vergence of the Navier-Stokes equation and using (|24|). Substituting (|27|) in 
the expression for the pressure in Proposition 1 one obtains the evolution 
equation 

\u\ 2 

Tn = RiRjiu 1 ^) -~Y + C ( 28 ) 

for n. Incompressibility can be enforced either by solving at each time the 
static equation (p5|) or by evolving n according to ([28]). 



Proposition 2. Let u be given by (^) and assume that the displacement 
solves ( jl^ and that the virtual velocity solves ([7^/ Assume in addition that 
the potential obeys j{2di ) (respectively (\23j)). Then u obeys the incompressible 
Navier-Stokes equations, 

d t u + u ■ Vw — uAu + Vp = /, V • u = 0, 

the pressure p satisfies $27 ) and the potential obeys also (\2B{) (respectively 
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The same results hold for the case of the whole R 3 with boundary condi- 
tions requiring u and £ to vanish at infinity. In the presence of boundaries, 
if the boundary conditions for u are homogeneous Dirichlet (u = 0) then the 
boundary conditions for v are Dirichlet, but not homogeneous. In that case 
one needs to solve either one of the equations fl25|) , (|28|) for n (with Dirichlet 
or other physical boundary condition) and the v equation ([Tj^) with 

v = V An 

at the boundary. 

Proposition 3. Let u be an arbitrary spatially periodic smooth function and 
assume that a displacement £ solves the equation ( |T3[ ) and a virtual velocity 
v obeys the equation ( \TQ ) with periodic boundary conditions and with C com- 
puted using A = x + £ . Then w defined by 

Wi = (<9^"> m (29) 

obeys the cotangent equation 

Tw + (Vu)*w = f. (30) 

Proof. The proof is a straightforward calculation. One uses ( |TTD to write 

Twi = (cU m )I\; m + v m T(d t A m ) - 2^0^)8^. 

The equation (p0|) is used for the first term and the equation ( |T5D for the 
second term. One obtains 

TWi = fi- (d iUj ) Wj + 2U {{d^Cr^d.Vr - (d k d t A m )d k V m } . 

The proof ends by showing that the term in braces vanishes because of the 
identity 

(diA m )C r ^ q]m = d q diA r . 



An approach to the Euler equations based entirely on a variable w (JTT 



([12]) is well-known. The function w has the same curl as u, u> = V x u = 
V x w. In the case of zero viscosity and no forcing, the local helicity w ■ u is 
conserved D t (w-uS) = 0; this is easily checked using the fact that the vorticity 
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obeys the "tangent" equation D t u = (Vu)u and the inviscid, unforced form 
of ([BOD. The same proof verifies the Kelvin circulation theorem 

(f w ■ dX = 

dt J 7 (t) 

on loops j(t) advected by the flow of u. Although obviously related, the two 
variables v and w have very different analytical merits. While the growth of 
w is difficult to control, in the inviscid case v does not grow at all, and in 
the viscous case its growth is determined by the magnitude of C which starts 
from zero. This is why we emphasize v as the primary variable and consider 
w a derived variable. 



6 Gauge Invariance 



Consider a scalar function (p. If one transforms v i— > v = v + Va4> an d 
n I— ► h = n + (f> then u remains unchanged in (JTJ) : uhu. The requirement 
that Ve • u = does not specify this arbitrary <fi. 

Assume now that the scalar cf> is advected passively by u and diffuses with 
diffusivity v. 

T(j) = 0. 

Then, in view of QI7D , if v solves fll9f) then 

v = v + V 'a4> 
also solves (|T^). If n solves then 

n = n + 

also solves fl28l) . If w solves the equation (p0|) then 

u> = + Vg0 

also solves ( jS0|) . The vector fields obtained by taking the Eulerian gradient of 
passive scalars are homogeneous solutions of (ppp. The vector fields obtained 
by taking the Eulerian-Lagrangian gradient of passive scalars, Va0 are ho- 



mogeneous solutions of (|19|). This can be used to show that if one chooses 
an initial datum for v that differs from uq by the gradient of an arbitrary 
function O there is no change in the evolution of u. 
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Proposition 4. Let each of two functions Vj, j = 1,2 solve the system 



with periodic boundary conditions for £j = Aj — x. Assume that the ini- 
tial data for Aj are the same, £j(x,0) = 0. Assume that each velocity is 
determined from its corresponding virtual velocity by the rule 



Assume, moreover, that at time t = the virtual velocities differ by a 
gradient 

Pfi = ~Pv 2 = u . 
Then, as long as one of the solutions Vj is smooth one has 

Ui(x, t) = u 2 (x, t), Ai(x, t) = A 2 (x, t) 

The same kind of result can be proved for ( [30]) using the Eulerian gauge 
invariance. 

7 K-bounds 

We are going to describe here bounds that are based solely on the kinetic 
energy balance in the Navier-Stokes equation ( ([|15j) and references therein). 
These are very important, as they are the only unconditional bounds that 
are known for arbitrary time intervals. We call them kinetic energy bounds 
or in short, K-bounds. We start with the most important, the energy balance 
itself. From the Navier-Stokes equation one obtains the bound 



T( Uj ,V) Vj = 2uC j Vv j + Q* j f 



with periodic boundary conditions, coupled with 



r(uj, v)Aj = o 



Uj = 



P ((VAjTvj) . 



1 



u(x, t)\ 2 dx + v 




1 



Vm(i, s)\ 2 dxds < Kq 



(31) 



with 



K = min{/c ; &i} 



(32) 
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where 



k = 2 J \u(x,t )\ 2 dx + 3(t-t ) J J \f(x,s)\ 2 dxds (33) 

and 

fci = J \u(x,t )\ 2 dx + - J j \A-^f(x,s)\ 2 dxds (34) 

Note that we have not normalized the volume of the domain. The prefactors 
are not optimal. The energy balance holds for all solutions of the Navier- 
Stokes equations. We took an arbitrary starting time to- The bound Kq is a 
nondecreasing function of t — to- We will use this fact tacitly below. In order 
to give a physical interpretation to this general bound it is useful to denote 
by 

e(s) = uL~ 3 J \Vu(x, s)\ 2 dx 
the volume average of the instantaneous energy dissipation rate, by 



E (t) = ^j3 J W(x,t)\ 2 dx 

the volume average of the kinetic energy; for any time dependent function 
g(s), we write 

1 A* 

(#(•))* = T T- / 9 i s ) ds 

for the time average. We also write 



F 2 = (l-* j \f(x,.)\ 2 dx^ 

G2 = ( L ' 3 j \^f(x,-)\ 2 dx 
and define the forcing length scale by 

G 2 



L 2 



f p2 ■ 
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Then ( PT|) implies 

2E(t) + {t- t ) (e(-)) t < iE{t ) + (t- t )F 2 mm j^; 3(t - *„)}• (35) 
After a long enough time 

t-t >-t, 
3v 

the kinetic energy grows at most linearly in time 

„ (it- t )L 2 f 
E(t)<2E(t )+F 2 r- ^-J. 

The long time for the average dissipation rate is bounded 

F 2 L 2 f 

lim sup (e(-))i < — e B- (36) 

These bounds are uniform in the size L of the period which we assume to 
be much larger than Lf. If the size of the period is allowed to enter the 
calculations then the kinetic energy is bounded by 



.,L 2 rF 2 { n Jj 2 fF 2 \ v(t-tn) 

E(t) < L 2 ^- + (E(t ) - L 2 ^ ) 



where 



L 2 F 2 = snpL- 3 j \(-A)^f(x,t)\ 2 dx. 
This means that for much longer times 

L 2 

t-to> — 

V 

the kinetic energy saturates to a value that depends on the large scale. But 
the bound ( |35| ) that is independent of L is always valid; it can be written in 
terms of 

B = 4£(t ) + (t- t )e B (37) 

as 

E(t) + (t- t ) (e) t < B. (38) 
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A useful K-bound is 

-t 



/ \\u(-,s)\\ L oo {dx) ds < (39) 

J t 

The constant has dimensions of length and depends on the initial kinetic 
energy, viscosity, body forces and time. The bound follows by interpolation 
from ([16[]) and is derived in Appendix A together with the formula 



Koo = C |^ + y^To) + t -jr-J* \\f(; s)\\Us} . (40) 

The displacement £ satisfies certain K-bounds that follow from the bounds 
above and (|I3|). We mention here 

M(-^)\\L^(dx) < / IK', s)\\L°°(dx)ds < iQo, (41) 

The inequality (|41|) follows from (|13|) by multiplying with £\£\ 2( - m ~ 1 \ integrat- 
ing, 

—4 / \((x,t)\ 2m dx + u f \ve(x,t)\ 2 \e(x,t)\ 2( - m - 1) dx+ 

2m at J J 

I \V\£(x,t)\ 2 \ 2 \£(x,t)\ 2{m - 2 '>dx+ 



+ / u(x, t) ■ £(x, t) \£(x, t) | 2(m - 1} dx < 0, (42) 



and then ignoring the viscous terms, using Holder's inequality in the last 
term, multiplying by m, taking the m-th root, integrating in time and then 
letting Tn — > oo. 

The case m = 1 gives 



Ay Wl , t)|2dx+ ,y |w(M) |^<^/ Wl , t) |^ 

and consequently, we obtain by integration from t = 



IJ WM < t VK , (43) 
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and then, using we deduce the inequality 



\V£(x, s)\ 2 dxds < 



K e 

2u 



(44) 



Now we multiply ( |13|) by — A£, integrate by parts, use Schwartz's inequality 
to write 



|/|W(M)l^ + ,/lMM)l^^/|V M (M)I^J/|V 4 (M)l^ 



-2 J Trace {(V£(x,t))(\/u(x,t))(V£(x 1 t))*}dx 
and then use the elementary inequality 



\V£(x,t)\ 4 dx) <C\\£(-, 



\L°° 



\A£(x,t)\ 2 dx 



in conjunction with the Holder inequality and fl4T|) to deduce 
d 



dt 



\V£(x,t)\ 2 dx + u / \A£(x,t)\ 2 dx < 



/|W(M)lW/|V M (M)l^ + ^/|V M (M)l^. 



We obtain, after integration and use of ([H], |4]) 



|W(M)| 2 dx + i/ 



\A£(x,s)\ 2 dxds < C 



Kot , K 2 K 



+ 



Recalling the bound ([J7|, |38|) on kinetic energy we have: 



(45) 



Theorem 1. Assume that the vector valued function £ obeys (IS) and as- 
sume that the velocity u(x, t) is a solution of the Navier-Stokes equations ( or, 
more generally, that it is a divergence- free periodic function that satisfies the 
bounds and (ffity)- Then £ satisfies the inequality together with 



-Sg J \£(x,t)\ 2 dx < (4E + te B )r 



(46) 
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and 



mxMjji + v fj mx,s)\^s< C (^ + ^). (48) 
In t/tese inequalities 



E = -^jj J \u(x,0)\ 2 dx, 
B = 4E + te B 



and €b is given in pq/ 

Let us consider the pair dispersion 



(5 t 2 ) = L 6 / / | a: - y| 2 dxay (49) 

7 J{(x,s/);|A(x,t)-A(y,t)|<5 } 

Using the triangle inequality (§) in (f46[ ) we obtain 

Theorem 2. Consider periodic solutions of the Navier-Stokes equation with 
large period L, and assume that the body forces have Lf finite. Then the pair 
dispersion obeys 

(b~t) < 3<5g + 24tE t 2 + 6e B t 3 . (50) 

Comment Use of the ODE ^ = u(X, t) requires information about the 
gradient VA and produces worse bounds. 



8 e-bounds 

This section is devoted to bounds on higher order derivatives of i. These 
bounds require assumptions. We are going to apply the Laplacian to (|IBD, 
multiply by At and integrate. We obtain 

~ J \A£(x,t)\ 2 dx + v J \\/A£(x,t)\ 2 dx = 

J d k u(x, t) ■ d k A£(x, t)dx + I (51) 
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where 



Now 



where 



I = J d k (u(x, t) ■ V£(x, t)) ■ d k A£(x, t)dx. 

/ = / ( a lM) .w(M).^,^ + // 

// = / M (,, t ).V W (M)).A W M) rf x 
and, integrating by parts 

// = -/a lM (x,t).V(^(M)).W(M) rf x 

and then again 

11 = J diu(x, t) ■ Vdid k £(x, t) ■ (d k £(x, t))dx 

and so 

1 = J diu i (x,t)d k £ j (x,t){d i d k + 6 ik A}d l £ j (x,t)dx 
Putting things together we get 

|/| < c\\ve(;t)\\ LX \\vu(-,t)\\ L2 \\vAe(;t)\\ L 2 

Thus 

~ J \A£(x,t)\ 2 dx + u J \VA£(x,t)\ 2 dx < 

° J |V M (x,t)| 2 dx + C||V£(-,t)|Uoc||V M (-,t)|| L2 ||VA£(-,t)|| L2 (52) 

Now we use an interpolation inequality that is valid for periodic functions 
with zero mean and implies that 

||W|Uoc<c||A£||2 2 ||VA£||2 2 
Using this inequality we obtain 

d_ 

di 



/|MM)l^ + ,/|VA <M |^< 
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£ J \Vu( Xl t)\ 2 dx + Cv-*\\Vu(-MlAW-MI* (53) 
Therefore we deduce 

/ |A^,t)| 2 ^<c^exp{^£ e 2 ( S )^} (54) 

where 

e(s) = vL~ 3 / \Vu(x,s)\ 2 dx (55) 



is the instantaneous energy dissipation. 



Proposition 5. If £ solves ( \13j) with periodic boundary conditions on a time 
interval t G [0,T] and if the integral 



T 

2 



e 2 (s)ds 



is finite, then 



\A£(x, t)\ 2 ^ + uj* J \VA£(x, s)\ 2 ^ < c^expl^J* e 2 (s)ds 
holds for all < t < T . 



9 Bounds for the virtual velocity 

We prove here the assertion that v does not grow too much as long as the 
norm of C is not too large. We recall that v solves (|20|) 

Tvi = 2vC mMi d k v m + Qjifj. (56) 

We multiply by Vi\v\ 2<yrn ~ is> and integrate: 

\v(x,t)\ 2m dx + v [ \Vv(x,t)\ 2 \v(x,t)\ 2{m - 1) dx+ 



1 d 
2m dt 



/ \V\v{x,t)\ 2 \ 2 \v(x,t)\ 2 ^dx 
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2v 



J C m , k ^x, t) (d k v m (x, t))vi(x, t) \v(x, i) | 2(m_1) cb+ 



We bound 
2v 



+ j Q Jl (x,t)f 3 (x,t)v i (x,t)\v(x,t)\ 2 ^dx. 
J C mMi (x, t) (d k v m (x, t))vi(x, t)\v(x, t)\ 2{ - m -^d. 



< 



\Vv{x,t)\ 2 \v{x,t)\^ m ~^dx + v / \C(x,t)\ 2 \v(x,t)\ 2m dx 



where 



and we bound 



\c(x,t)\ 2 =J2\c m ,kA^t)\ 2 , 

m,k,i 



Qji(x, t)fj(x, t) Vi (x, t) \v(x, t)\ 2(m -^dx 



< 



\g(x,t)\ 2m dx} <j / \v(x,t)\ 2m dx 



1 2m, 



where 



9i(x,t) = Qji(x,t)fj(x,t). 
The inequality obtained is 

<2mu J \C(x,t)\ 2 \v(x,t)\ 2m dx + 



+ 2m( I \g(x,t)\ 2m dx} <{ / \v(x, t)\ 2m dx 



1 2m, 
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Let us consider for any m > 1 the quantity 

q(x,t) = \v(x,t)\ m . 

The inequality ( j60| ) implies 

lj(q( x ,t)f + 4u(l-^-) | |Vg(x,t)| 2 <2mz, J \C(x,t)\ 2 (q(x,t)) 2 dx 



+ 



+ 2m|y |^(x,t)| 2m c/x| 2m |y (g(x,t)) 2 dx 

Using the well-known Morrey-Sobolev inequality 

i 

(q(x)) 6 dx^ <C { I \Vq(x)\ 2 dx + L~ 2 I (q(x)) 2 dx 



and Holder's inequality we deduce 



j t [( q (x t t))* + 4v(l-±) I |Vg(M)| 2 < 



<2muC { I \C{x,t)\*dx\ <j I \Vq(x,t)\ 2 dx + L~ 2 I (q(x,t)) 2 dx 



+ 2m jy t) \ 2m dx | | y (g(ar, t)) 2 dx 

Therefore, if on the time interval t G [0, r] C(x,t) obeys the smallness con- 
dition 



\C(x,t)fdxY< J 2 -^J± (61) 



then we have the inequality 



j t H-,t)\\ L 2m < U ^ 2l ^ \\v(-,t)\\L^ + \\g(-,t)\\L^ 



for t G [0, r] and consequently 



i/(m-l)t /"* 

K-,£)||£2m < ||v ||Lame ^ + / \\g(-,s)\\ L 2m (62) 







holds on the same time interval. 
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10 Appendix A 

In this appendix we prove the inequality (^) and derive the explicit ex- 
pression for Koo. The calculation is based on ([0). All constants C are 
non-dimensional and may change from line to line. Solutions u of the Navier- 
Stokes equations obey the differential inequality 

-7- f I Vw(x, s)\ 2 dx + v [ I Au(x, s)\ 2 dx 
ds J J 

< ^ (J \Vu{x,s)\ 2 dx^J + ^- J \f{x,s)\ 2 dx. 

The idea of (Jl6|) was to divide by an appropriate quantity to make use of 
the balance (pTj) . The quantity is 



(G(s)) 2 = (V + f \Vu(x,s)\ 2 dx^j 



where 7 is a positive constant that does not depend on s and will be specified 
later. Dividing by (G(s)) 2 , integrating in time from to to t and using (|3T| ) 
one obtains 

'* \Au(;s)\\UG(s))- 2 ds< 



is 5 z/7 2 i/ 2 7 4 J to 



The three dimensional Sobolev embedding-interpolation inequality for peri- 
odic mean-zero functions 

1 1 

|| u || ^oo < C|| Vm||? 2 II Am||? 2 
is elementary. From it we deduce 

||«MIU- < C\\Vu(;s)\\l 2 (G(s)) 1 * [HA^-,,)^^^)- 1 ]^ 

Integrating in time, using the Holder inequality, the inequality (|3~1~D and the 
inequalities above we deduce 

t 

\\u{-, s)\\ L oo < Cr 

to 
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where the length r = r(7, t, u, K ) is given in terms of six length scales 

K v 2 (t - 1 )7 2 

—f = r , -= = n, = r 2 , 

z/ z 7 Z v 



to 



v 

and 

»"5 = vK^ - *o)- 

The expression for r is 



r = r + (r )4(n)4 + (r )2(r 2 )2 + (r ) 4(^)4 + 



111 31/ 7*1 \ 4 

(r )3(r 4 )3(r 5 )2 + (r )4(r 4 )4 



>'2 



The choice 

7 4 



z/ 3 



t-* 
entrains 

ri = r 2 = r 3 = r 5 

reducing thus the number of length scales to three, the energy viscous length 
scale r , the diffusive length scale r 5 and the force length scale r 4 . The bound 
becomes 

Koo = C(r + r 4 + r 5 ) 

i.e. fHCh. 



11 Appendix B 

We prove her the commutation relation ([17]). We take an arbitrary function 
g and compute [r, L^g]) where Y = T u (u, V) and Li = V^. We use first (flTD : 

[r, Lfflr] = r (Qjidjg) - QjidjTg = 

T(Q ji )d j g + QjiTdjg - 2vd k (Q ji )d k d j g - QjidjTg = 
(commuting in the last term dj and T) 

F(Qji)djg - 2vd k {Qji)d k djg - Qjidj{u k )d k g = 
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(changing names of dummy indices in the last term) 

(T(Qji) - Qkid k {uj)) djg - 2vd k (Q ji )d k d j g = 

(using (H)) 

2uQ jp (did k A p )(d k Qii)d j g - 2ud k (Q ji )d k d j g = 

(using the definition (||) of Va 

2u(did k A p )(d k Q H )(L p g) - 2ud k (Q ji )d k d j g = 

(renaming dummy indices in the last expression) 

2v(d k (Qn)) {(did k A p )(L p g) - did k g} = 

(using (0) in the last expression) 

2v(d k (Q u )) {(did k A p ) (L p g) - d k (di(A p )(L p g))} = 

(carrying out the differentiation in the last term and cancelling) 

^(OkiQuMiA^dkiLpig)) = 

(using the differential consequence of the fact that Q and VA are inverses of 
each other) 

2vQ u {d k di{A p ))d k {L p g) = 
(using the definition (H) of V^) 

2uLi(d k (A p ))d k (L p g) = 

(using the definition (||) of C mtk]i ) 

2vC Ptk .id k (L p g) , 

and that concludes the proof. We proceed now to prove fll8|). We start with 

(ED 

T(d k A m ) = -(dhUj^djAm) 

and apply Lf 

L t (T(d k A m )) = -Li {(4%)(^A m )} . 
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Using the commutation relation ( p]) and the definition we get 

r(C m , fc ;i) = -Li {(dkU^idjArn)} + 2vC p ^ i d l L p (d k A m ). 

Using the fact that Lj is a derivation in the first term and the definition (|9]) 
in the last term we conclude that 

r(Cm,fc;i) = —(dkUj)C m j-i — (8j A m ) (Lj(9fcMj)) + 2vC P) l-i{dlC m fop) 

which is (0). We compute now the formal adjoint of 

{V\Tg = -djiQiig) = -V\{g) - {d^Q^g 

(with ©) 

(y^Yg = -V l A (g) - {(d 3 A p )L p (Q Jl )}g = 
(using the fact that Q is the inverse of VA) 

{V\y g = -V\(g) + QaC^g. 
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